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I— I . Abstract 
^ ■ 

CIh' DGLAP evolution equations are modified in order to use all the quark families in the full 

scale range, satisfying kinematical constraints and sumrules, thus having complete continuity 
D I for the pdfs and observables. Some consequences of this new approach are shown. 

Comments: 12 Pages and 5 Figures 

> ■ 

^ 1 Introduction 
(N 

• As it is well known heavy quarks present a challenge in the phenomological description of deep inelas- 
tic scattering of leptons against nucleons: DGLAP [1] evolution equations are essential ingredients 
to this description but they consider only massless partons. With the advent of HERA results, the 
increase of the statistical precision of the measurements and the prespectives open by LHC, the 
necessity to have a way to predict the observables from the lowest to the highest scale became 
^ ■ clear and several so-called variable flavor number schemes (VFNS) have appeared: 

• In the massless approach a flavor h is usually considered to be active when its quark mass 
is smaller than the scale . Up to this limit its parton distributions (for quark and antiquark) 
are null, after the number of flavors is increased by one and they start to evolve from zero. 

• In the massive approach the number of flavors is considered to be constant but the heavy quark 
is produced by the existing partons thus it have structure functions coming from dependent 
coefficient functions and no (pdf) distributions. 

• VFNS approaches make use usually of the two preceeding approaches according to predefined 
regions: At some arbitrary fixed value the number of quark fiavors is increased by one and 

the pdfs are rearranged according to theory inspired rules before resuming evolution. All this 
is a still an ongoing activity. 

The path of this work is different: It will consider always six quark species but modify coherently 
DGLAP equations, splitting functions, coefficient functions,fiavor number and running in order 
to have all the parton distributions and structure functions continous in the full kinematical range. 
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This will be done furthermore in such a way that when the phase space for an heavy quark is quasi 
null {ml ^ ) its distributions are quasi null and its influence on the full system is also quasi null. 
On the contrary when its phase space is quasi complete (m^ <C ) it will behave like a light quark. 
The key to this will be kinematic constraints and sumrules. 

To fulfill these commitments, section 2 presents DGLAP equations in an appropriate manner, enpha- 
sizing the components which will need modification. Section 3 will present the modified (extended) 
equations defining the method thereafter called cfns (continous flavor number scheme), section 4 
brings arguments to justify the method, section 5 shows a possible program implementation and 
finally Section 6 will present a comparaison between massless scheme and cfns at NNLO in a^. 



2 Kinematics and DGLAP equations 

The electron proton reaction is: 

e{l) + P{p) e{l -q) + X{p + q) 

Within the symbols ( ) are quadrimomenta. 

Kinematic variables arc: 

S=il + p)' Q' = Sxy = -q'>0 x = g y = g ^ {p + qf = + Q^'- - 1) 

M and W are the initial and final hadronic masses. 

Parton o may be kicked out of the target if the final hadronic mass is W > 2mo where rrio is the 

parton o mass. This translate into a kinematic limit x <lo with lo—{l-\- A.m'^Q~'^^ 

Light partons fulfill always this condition but heavy quarks do only for Q"^ ^ oo. 

Usually the flavor number A/} is taken as the number of quark families such that > rri^ : A/} = Nf 

integer. A// is the main concern of the overall approach. 

DGLAP equations read : 

V are splitting functions, i and o are parton distributions (pdf) and run on the 1 + 2A^/ partons 
species, they will be noted also by the name of their species p — g,d,d,u,u,... and for the quarks 

= d±d,... will be introduced. 
All these are x and functions, (g) note the convolution between two functions of x defined by: 



[A ® B]{x) = £ A{z)B{^J^ = £ Ai^)Bit)^ 



The following properties hold: 

A^B^ B^A [A^ B]^C ^ A<^[B ^C] O S] = [xM] ® [a;"S] 

J^[A ^B]dx^ Jq A{x)dx B{x)dx 

/o 5{x - l)dx = 1 x5{x - l)dx 5{x-l)®A^ \A{^) 
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DGLAP equations separate into two independent subsystems the first one being: 

1^ = V,, ® gm + ® ^HQ') (2) 

q=d 



Using definitions inspired by [2]: V%s = K ± ^J?' ^5 = ± ^fr = jfj^^t = J^- 

The second system may be deduced from the first by replacing superscript + by — and suppressing 

all references to gluon g. So it will not be mentioned anymore and even superscript + will not be 



written in the following 
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Notice that V^g existence comes from the /3o term in Vgg. 
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3 Modified DGLAP equations 

The idea is to modify the kernels, keeping them functions of a single argument (apart from de- 
pendance coming from A// and a^) in order to satisfy simultanously the three kinematical constraints 

Voi gives the change to outcoming parton o at Bjorken x radiated by incoming parton i at Bjorken 

X 

z ' 

Problematic cases are when parton o is heavier than parton i like for c ^ 6 as shown Figure [1] where 
the corresponding term in equation(IT]) is depicted. For a given value of x the convolution integral 
variable z runs vertically inside the big triangle at least for a massless parton. But for a massive 
quark the rightmost triangle has to be removed and if one wants to keep the splitting nature of the 
convolution only the leftmost hached triangle has to be kept. Notice that this will also supress an 
unwanted discontinuity of that term aX x = lo and bring it gently to 0. 

An other way to present the modification is to replace P by 5 in the problematic changing term: 

jyo.ijW)j (4) 

Requesting this term to be null for x > Iq means that S{u) = for u > lo which is satisfied by: 

Soi = Voi ® ICoi (5) 

With the definition JCoi = foiS{x — Ig), where foiiQ"^) is not decreasing and goes to 1 for -> oo 
Notice that the effect of the 6 function is to replace V{x) by P(0 with x = E,lo 
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With this modification the subsystem becomes: 



1^ = V,, ® K,,, ® gm + J] P,, ® K,,, ® qm (6) 



q=d 
t 



r=d 



3.1 Sumrule constraints 

In the above equations P's are functions of the A// to be defined now. The momentum sumrule 
imposes constraints which have to be satisfied (and are for the standard DGLAP). They are obtained 
by requesting that the first momentum of the sum of all parton distributions is constant and equal 
to 1 for any value. Introducing A = JCxdx and Q = Vxdx, taking the derivative of this 
momentum, using the DGLAP equations and the property that the n-momentum of a convolution 
product is equal to the product of the n-momenta of its components, one find easily: 

2..A.. + ^EN. = (8) 

^ q=d 

Q,q\q + QnS^Z' + ^ E ^'■^ = (9) 

^ r=d 

Using the fact that those equations are satisfied by the original DGLAP equations where all the As 
are 1, one gets: 

.r _ T.q=d\9 _ Qgq + QnS . „s 

^gg ^gq^gq + ^NS^qq ^ 

Leading to the final extended DGLAP equations Q 

= V,, ® gix, Q') + j^V,,® q^x, Q') (11) 

q=d 

= ® [K ® 9ix, Q') + Y.VI® r+(x, Q')] + ® qH^, Q') (12) 

T=d 

Where K,q = 6{x — Iq) . Furthermore Q{V) being polynomials in the flavor number, the only solution 
is to redefine the latter as being Aff = Yll=d h- 

Quark kinematical limit or longitudinal phase space Ig then may be also viewed as quark activity 
going smoothly from to 1. 



^ The simplest hypothesis has been used: fio — Agg = A^'^ = 1 
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3.2 Subsystem decoupling 



As it is well known the usual DGLAP equations decouple in a g,!! system of two coupled equations 
and several non-singlet independent equations. It is now not so simple. The g,Yl singlet system 
still decouple but the non-singlets cannot all decouple. In fact there is not anymore a unique way 
to simplify the system. The following solution is one of the both which have been realised in the 
programn implementation section 5. 
Define: 

Subscript ln is used do distinguish these non singlets from the usual ones which may not be used 
here with A// varying continously. With these one get the following subsystem: 



2^ 



dg{Q 
dln{Q^] 

dln{Q^] 

dlLNm 
dln{Q^) 



= Vgg(^g{Q^)+Vgg^^Q^) (13) 

= /C (2) (8) ^(Q') + + /C (g) 7^5) <H) S(Q2) (14) 

= Vns®Iln{Q^) (15) 

= }Ch®Pgg^g{Q^)+rNs®h{Q^) + ]Ch^Vs®^Q'') (16) 



dln{Q^) 

After evolution of the seven pdfs g,I^,c,b,t,dLN,ULN the full system may be recovered using: 
Ei = E- c- 6- t dLN + uln + sln = 

Note that for //j — >■ the corresponding DGLAP equation will get decoupled (Appelquist-Carazzone 
theorem) and the kinematical constraint automatically verified. 



3.3 Renormalisation equation 

As seen above the momentum integral sumrule leads to a specific non integer value of A// and as a 
consequence also for /3o and by extension for the full set of f3 governing the ag running. It is also 
natural that the coupling constant depends on the sum of flavor activities Ig and not only on flavor 
number. 

3.4 Coefficient functions 

As it is the structure functions and not the parton distributions which are observable one has to find 
also a procedure to modify the coefficient functions. 

• The transform parton distribution — > structure function has exactly the same structure that 
the one of DGLAP equations, it is obtained by the changes: 

do{Q') ^ ^ 
dln{Q^) 
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• This transform is in fact nothing more than a change of scheme, an example is going from MS 
to DIS for F2. 

For Fi and F3 the schemes are unnamed but they still exist. 

From this one may infer that coefficient functions have to be modified in the same way that 
splitting functions. 

In fact it is the importance of kinematic constraints stressed by one of the R.Thorne papers and its 
use in the latest schemes which started this work. 
Charged currents 

Exactly the same procedure will be used, the phase space only will change using: 

/o = (l + m2Q-2)-i (17) 

4 Theoretical Considerations 

Extended DGLAP equations have been set up very close to the ordinary ones. 

• They reduce to it when scale is far away from any heavy quark mass squared. 

• They are integro differential equations, linear in all the pdfs which they use in conjonction with 
convolution integrals. 

• They satisfy the kinematic constraints of the heavy partons. 

• The resulting pdf are continous in x and at least if the input pdfs at the initial scale are. 

• The derivatives of the pdfs are also continous. 

• The same properties apply for charged and neutral currents structure functions Fi,F2 and F3 
(see Figure HI). 

• When there is only active fiavors (all Ih close to or 1) Soi{x, Q"^) Poi{,x, Q'^,J\ff = Nf) 

• They are supported by works and concepts which are not really new: 

— ^ the scaling variable was used in many papers, see [3] as an example. 

— [3] use also anomalous dimensions variable with Q"^. Anomalous dimensions leading 
to splitting functions their arguments should hold here. They advocat^ li ~ (1 + 

— Moreover |3] presented a /3-function variable with (due to /3q ). They advocate here 
h ~ (l + 5m2Q-2)-i. 

— [1] have considered A// to be and even order dependent and there is still develope- 
ments on the /3-function [5]. 

— Small are the differences between massless and continous behavior for (see Figure [2]). 

"^they did note have the same name and used it differently however. 
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— The procedure leading to satisfaction of the kinematical constraints as been used lately 
for coefficient functions in GM-VFNS schemes |7j. 

— Presented here is the simplest solution to the chosen goal but it is still possible to make 
modifications making use of the freedom brought by the /oi(Q^)s following |1] example. 

Satisfiying all the points it was designed for, it as some serious advantages: 

• It does not mix up different as orders as do mixed schemes. 

• Heavy quarks participate to the evolution when they start to appear, that is at the beginning 
and at very low x and even at leading order in as- 

• There is only internal partons, no externals. 

• It covers the charm-bottom region where they are both opening up (see A// graph Figure [2]) 
which is not yet the case in VFNS. 

• It should give a better treatement of the small x region where heavy quark pdfs appear first. 

It seems that ressummation is done, due to evolution, but not the forward divergence pole subtraction. 
In that case this method should be considered only as a phenomenological model which may be 
optimized by adjusting the foi functions with the help of the rates given by [8j. The comparaison 
between the F^'^ of the two approaches would however be not that simple, having so different concepts 
(quark internal versus external). 

An important point is that nature has 6 flavors and not 3,4,5 or 6 depending on physicist will, so 
neglecting that fact is making an approximation. This paper is also making approximations which 
are certainly valid when A//(Q^) is close to an integer. Maybe renormalisation group theory could 
bring some light to this. 

5 Program implementation 

Modiflcations to be made to existing programs may be done without changing their conceptual build- 
ing, however changes are not that small. The QCDFIT case ([6]) will be presented in some details 
as an example. 
QCDFIT concept 

It is a program which works in x space. It includes an optimizing procedure (Minuit package [TO]), 
it accept a variety of input distributions. It has a variety of outputs: Pdfs, cross section for lepto 
production ,Drell-Yang mecanism... It precalculate the full evolution saving there a lot of CPU time. 
Pdf representation 

QCDFIT uses two grids: The x grid has linear spacing in log(x). There is much more liberty for the 
grid. Present version uses a grid approximatly equally spaced in a^. [fl All the functions of 
X or or (x, Q"^) are calculated on these grids. To know values elsewhere a linear approximation 
has to be made (there is specialised routines for that). One consequence is that distributions are 
represented by vectors and that the evolution of the distributions from one grid value to an other 

maximum of 20000 x nodes and 600 nodes has been used 
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q2 gj,j^ value is given by matrices. 
Kernel representation 

They are linear operators applied to the pdfs and so after discretisation have a matrix representa- 
tion. But due to the linear appoximation used, the structure of convolution and more deeply to the 
concept of parton branching those matrices are upper triangular band matrices: their elements are: 
J^ij = rrii-j with i > j. 

As a consequence in QCDFIT they are represented by one dimensional arrays and a system of fast 
and simple routines have been developed to deal with kernel multiplication, inversion, exponentia- 
tion, square root and vector (pdf) multiplication. Notice that target mass correction, higher twist 
and renormalon a la M.Dasgupta and B.R. Webber are eventually computed in this kernel frame. 
QCDFIT Evolution 

Integration of the renormalisation group equation is made analytically for the needed flavors. Inte- 
gration of a" needed for kernel integration is made at the same time and with the same technique. 
For the NS case it is easy to show that the solution of DGLAP involves an integration of the kernel 
followed by an exponentiation (easy and fast operation with QCDFIT concept). For the Singlet case 
things are more complicated and an additional calculation is needed. 
QCDFIT Modification 

The representation of the Pdfs, the kernels, the convolution algebra stay the same, but the evolution 
itself has to be modified: Integration of the renormalisation group equation is now made numer- 
ically as its /? parameters are functions of A// and so of Q^. The transport matrices defined by 
o{Q'j_^_i) = '^iT{o,i) ® i{Q'j) are obtained by integration of the subsystems using: 



In the product Sln{Q^) has to be small enough to see only the rounding errors when increasing the 
number of nodeso 

6 cfns massless comparaison at NNLO 

The aim of the exercise is to show what kind of new features might be seen on cfns and how far they 
extend away from the transition points of the other schemes. 
Kinematical range used 

Very often the start of evolution Qlnput is chosen just below the charm masss squared in order to define 
pdf inputs only for light partons. But for cfns heavy quarks are always present if the kinematical 
range permits, so in order to have only light quarks present , a Qj^put — 0-6 Gev^ is used, low enough 
to justify neglecting all the heavy quarks at input. Needless to say that at so low a the predictive 
power of the pdfs is completly absent but it is a parameterless way to get a sensible charm when out 
of the non perturbative region. 

A data sample made of about 1850 F2 or cross section measurements extracted from NMC, BCDMS 
on protons and HI preliminary is used to fit the input pdfs independently for cfns and for massless 

'^QCDFIT use presently grids from 10000 to 500000 nodes 




(18) 
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scheme. The fitted distributions are g,Uvai, d^ahU = d = Xs with x dx s = 0.53 x dx d at 

g2 = lOGev^ [9] 

Results 

Figure [3] shows the momentum integral fractions for the partons species in both schemes as a function 
of log(Q^). Note for heavy quarks the earlier start and slower rise of cfns. Figures H] shows the cfns 
and massless structure function J-'g at x = 0.01 versus Q^. The three kinks at in massless 

are due to the alternative coming in of quarks and antiquarks wih negative sign. In cfns slope 
continuity is restored because there in a single DGLAP system and coefficient functions with no 
discontinuities. Figures [5] shows the cfns and massless top distributions at = 39800 where both 
momentum fractions are equal, cfns has a much steeper distribution with highers values at small 
X coming from its early developement and from the kinematical constraint (not shown is the same 
effect for charm and beauty). 

Conclusion 

Effects shown are impressive and justify the idea of evo luting heavy quarks even when they have to 
be considered as heavy. 

Work will continue to assess on more firmer theoretical ground these ideas. 
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Figure 1: C ^ B a problematic changing term. Convolution domain of Equation HI Vertical thin 
spaced lines show the kinematically reduced domain. 
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Figure 2: scheme comparaison with as{M'^) = .118 a) A/} versus \ogQ^ for cfns and 

massless. b) Relative difference between massless and cfns for versus logQ^- 




Figure 3: From top to bottom g , Umi , d^ai , u = d,s,c,b,t momentum fractions versus logQ^. 
cfns: black full, massless: red dotted. 
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Figure 4: T\ versus log cfns: black full, massless: red dotted. 




